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ABSTRACT 

Stars orbiting within 1" of the supermassive black hole in the Galactic Centre, 
Sgr A*, are notoriously difficult to detect due to obscuration by gas and dust. We 
show that some stars orbiting this region may be detectable via synchrotron emission. 
In such instances, a bow shock forms around the star and accelerates the electrons. 
We calculate that around the 10 GHz band (radio) and at 10^“^ Hz (infrared) the 
luminosity of a star orbiting the black hole is comparable to the luminosity of Sgr A*. 
The strength of the synchrotron emission depends on a number of factors including 
the star’s orbital velocity. Thus, the ideal time to observe the synchrotron flux is when 
the star is at pericenter. The star S2 will be ^ 0.015" from Sgr A* in 2018, and is an 
excellent target to test our predictions. 
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1 INTRODUCTION 

Over 100 young massive stars inhabit the central parsec 
of the Milky Way (for a review see iGenzel et al.l l2010l : 
iMapelli &: Gualandrij l2015l) . The stars whose orbit lies 
within 0.04 pc fro m the Galactic Gen t re (GC) are know n 
as the S-stars (e.g. ISchodel et al.l [ioo^ : iGhez et '^! 20051). 
The o rbits of 28 S-stars were determined bv iGillessen et al.l 
bOOgI) . 19 members have semimajor axis a ^ 1 arcsec. Of 
those, 16 are B stars and the rest late-type stars. Of partic¬ 
ular interest is the star S2 (also known as SO-2 ) which has 
been observed for more than one complete or bit ( Ghez et al.l 
I2OO5I : iGhez et al.ll2008l : iGillessen et al.ll200^ ). S2 orbits the 
supermassive black hole, Sgr A* * * § , every 15.9 years. It is 
a BO-2.5 V main sequence sta r with an estimated mass 
of -^15 Mq (iMartins et al.ll2008l ). The second complete or- 
bit of a star aroun d Sgr A* was announced not long ago 
llMever et al.l 1201^ 1. This star, S102 (also know as SO-102) 
has a period of 11.5 years and is about 16 times fainter than 
S2. Both S2 and S102 provide compelling evidence that Sgr 
A* has a mass of 4 x lO^M©. 

The detection of young, m assive stars with o rbits close 
to Sgr A* was surprising (eg. iGhez et al. |[20 o 3). A num- 
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ber of possible mechanisms have been propo s ed to account 
for the origin of the S-stars. iLockmann et ah] bOOSi ) argued 
that the dynamical interaction of two stellar disks in the 
centr a l pars ec could lead to the formation of the S-stars. 
Ig^ proposed that the S- stars were born in the 

disk a nd migrated inward. Recently. IChen fc Amaro-Seoanj 
b014h theorized that a few Myr ago the disk extended down 
to the innermost region around Sgr A* and Kozai-Lidov- 
like resonance resulted in the S-stars. Perhaps the simplest 
and arguably most likely scenario is that at least some of 
the S-stars are the result of a three-body interaction with 
Sgr A* (e.g.lGins burg fc Loebll200 d:lGinsburg fc Loebll2007l : 


_ uy ; _ 

iGinsburg et al. 20121 : Zhang et al. 20131 ). In this scenario, a 
binary star system interacts with Sgr A*, and tidal disrup¬ 
tion leads to one star falling into the gravitational well of the 
black hole w hile the com panion is ejected as a hypervelocity 
star (HVS) ( Hill j 19881 ). HVSs were first observed in 2005 
dBrown et alj [200f ) and as o f today som e 24 have been iden - 
tified (see Brown et al.ll20l3 for the list). lBrown et all boisl ) 
studied 12 confirmed HVSs and found that the vast majority 
are consistent with having a GC origin. Furthermore, obser¬ 
vations indicate t hat these HVSs are l ikely massive slowly 
pulsating B stars (IGinsburg et al.ll2013l) and thus consistent 
with the known S-stars. 

There are various candidate HVSs which a r e far 
less massive than B-type star s ( Palladino et al.l l2014l : 
IZhong et al.l l2014l : IZhAg et al.l l2015l) . Given ^e fact 
that the stellar density in the central parsec is ~ 
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10®Mq pc ^ dSchodel et alJlioogf) a distribution of masses is 
expected. However, detecting stars at the GC is notorio usly 
difficult due to dust extinction fe.g. ISchodel et ^l 201 flh . In 
this paper we discuss using synchrotron radio emission to 
observe and possibly monitor stars at the GC. Of particular 
importance is the fact that S2 will reach periapse in 2017. 
We show that the closer S2 is to Sgr A* the more likely we 
are detect the synchrotron emission. Thus, S2 serves as an 
ideal test subject. In Section 2 we provide the physics behind 
the synchrotron emission. In Section 3 we show that a star 
such as S2 may be observable via its synchrotron emission. 
We conclude with some discussions and future observations 
in Section 4. 


2 SHOCKS AND SYNCHROTRON EMISSION 


We consider stars orbiting Sgr A*. Interactions between 
winds from such a star and the interstellar medium (ISM) 
will create two shocks, a reverse and forward shock. The 
forward shock propagates into the ambient medium which 
is swept up and subsequently compressed and accelerated. 
The faster wind is compressed and decelerated by the re¬ 
verse shock. We are interested in the emission from electrons 
accelerated by the forward shock. In our case the forward 
shock is a bow shock with Mach angle 6 ^ , where A4 

is the Mach number. In such a scenario, roughly half the 
star’s mass loss contributes to the shock. The mechanical 
luminosity is simply given by kinetic energy which depends 
upon the mass loss rate, M, and wind speed 

Lw ~ (1) 

For a massive star such as S2, typical values for Mu, are 
~ 10~® Mq yr“^, Vw ~ 1000 km s“^ and thus is around 
10^® erg s~^. Consequently, the total non-thermal luminos¬ 
ity, Lnt, is given by the fraction of electrons, e„t, accelerated 
to produce non-thermal radiation 

Lnt — ^ntL w • ( 2 ) 

We let tnt be 5% although sim ulations show that this value 
is uncertain JGuo et al.l 2014l ~). The energy density of the 
amplified magnetic field is given hy Ub = /Stt. Therefore, 
assuming equipartition of energy, Ub = ^BnskTs leads to 

B = {Sn^BUskTs)^^^ (3) 


where ^b is the fraction of thermal energy in the magnetic 
field, Ua is the post-shock number density, and kTa the tem¬ 
perature of the post-shock medium. In the strong shock limit 
for an adiabatic index of 5/3, the post-shock number density 
ris is 4 times the number density of the ambient ISM. Al¬ 
though the value of is hi ghly uncertain, a value of 0.1 is 
reasonable llVolk et al.ll2005l l. The resulting values for B are 
~ 10“^ — 10“® G. The Rankine-Hugoniot jump conditions 
give the post shock temperature 

" Gt + lYM^ ° ^ ^ 

where To is the upstream temperature, and qt is the adia¬ 
batic index which is taken to be 5/3. At the contact discon¬ 
tinuity, ram pressure from the star/ISM is balanced by ram 
pressure from the wind. Thus we have p»u* = pwv'^ where 


the mass flux is given as p-wVw = 4 ^^ this leads to the 
standoff radius 


Ro 


, IT'ujV'u, 

M7rp*u* 


n1 /2^ 


(5) 


p* = noirip where nip is the proton mass. At a distance of 
around 0.3 pc from Sgr A* Uo ~ 10® cm~® an d the tem¬ 
perature of the ISM is ~ 10^ K llQuataertll2004h . Assuming 
V, = 1000 km s“^ a typical standoff radius is approximately 
10 “® pc. 

We consider a broken power law distribution of elec¬ 
trons generated via Fermi acceleration, written as 

NOd-y = Ko7“^(1 -I- ^ 7 < 'fmax) (6) 

"lb 


where Ko is the normalization factor in electron density 
distribution, p is the electron power law distribution index 
which in our calculations is ~ 2. 7 is the Lorentz factor, 
7 miTi and y-max are the minimum and maximum Lorentz 
factor respectively, and 75 is the break Lorentz factor due 
to synchrotron cooling. The total power from synchrotron 
emission of a single electron is given by 

P = LI^Vb^ (7) 

where ro is the classical electron radius and /3 = v/c 
liRvbicki fc LightmarJl98^ . The corresponding synchrotron 
cooling time is 


icool 


-ymc? 

P 


(8) 


The break Lorentz factor due to synchrotron cooling is ob¬ 
tained by equating the synchrotron cooling times scale and 
the dynamical timescale 



'^W 


(9) 


where Ro is the standoff radius and u™ the wind veloc¬ 
ity. The peak luminosity lasts for the pericenter crossing 
timescale with typical values ~ 1 yr. The precise shape of 
the peak depends on the orbital parameters, the ambient 
gas distribution, and the wind mass loss rate. The charac¬ 
teristic timescale for the peak is the pericenter crossing time, 
~ h/v. ')min is Set to One in our calculations, ymax is obtained 
by equating the acceleratio n timescale, tacc = GccRlc/v^ 
llBlandford fc Eichlerl Il987fl . to the dynamical or cooling 
timescale where ^acc is a dimensionless parameter on the or¬ 
der of unity. The Larmor radius is given by Rl = 'ynioC? jeB 
where is the mass of the electron and e is the electron 
charge. We find that taco ~ 10““^ yr, which is far shorter than 
the dynamic time. At radio frequencies the cooling time is 
a few orders of magnitude greater than the dynamical time 
and therefore synchrotron cooling is negligible. Thus, we as¬ 
sume that all the shocked electrons contribute to the ob¬ 
served synchrotron emission. 

The synchrotron flux and power are computed using 
the radiative transfer equation lIRvbicki fc LightmarJll986fl 
which leads to 


= ( 10 ) 

Olv 

where jv is the emission coefficient, Uq the absorption co ef- 
ficient, and the optical depth (see lWang fc Loeal2014l for 
further details). 

























Detecting Stars at the Galactic Centre via Synchrotron Emission 3 


year 



Figure 2. Synchrotron power (arcseconds) versus distance (erg 
s“^ cm~^) from Sgr A* for star S2. We let M = 10“® Mq yr“^. 
At a frequency of around 1 GHz, S2 is quite bright with a lumi¬ 
nosity of ~ fO mjy. The synchrotron power will be greatest in 
2017, when S2 is at periapse. 


3 EMISSION AROUND SGR A* 


Figure [T] shows the expected synchrotron flux for a star or¬ 
biting Sgr A*. For the left panel we kept the wind velocity 
constant at 1000 kms“^. We extrapol ated the particle den¬ 
sity to have the value 10“^ cm“® from lQuataeriJ ll2004l '). The 
mass loss ra te for hot massive stars is poorly constrained (for 
a review see lPuIs et aljfioo^ . Thus, we varied the mass loss 
rate between 10“^Mq yr“^ and 10~®Mq yr“^ which are 
acceptable va lues for hot massive stars such as the S-stars 
(lDupredl2015h . Given that 




vLv 1 
u 


( 11 ) 


where d = 8 kpc we get a fl ux between ^ 10 — 1000 mJy 
in the GHz range. Recently, lYusef-Zadeh et ^ ()2015l 'l de¬ 
scribed radio observations of over 40 massive stars within 
30” of the GC. Their values for the flux are consistent with 
our results. Similar results were obtained when we kept the 
mass loss rate constant at 1O~®M0 yr~^ and varied the wind 
velocity between 1000 km s“^ and 4000 km s“^ (see right 
panel of Figure [T]). 

In Figure [2] we plot the synchrotron power versus dis¬ 
tance from Sgr A* for S2. At periapse S2 has a speed of 
~ 5500 km s“^ while at apoapse the speed drops down to 
~ 1300 km s“^. A higher speed leads to a stronger shock. 
However, even at apoapse it may very well be possible to 
observe the synchrotron emission from S2. At 10 GHz, S2 
has a flux density of 10 mJy. Thus, S2 is an excellent 
target to observe. It is important to note that radio wave¬ 
length photons are scattered and the i mage size follow s a 
dependence (e.g. iBower et al.l l2006l : iFish et ^ l2014h of 
~ 1 mas (A /cm)^. Therefore, in order to resolve S2 at the 
0.01” level, we need to observe at a frequency of 10 GHz or 
greater. Sgr A* has bee n observed in the radio for decades. 
iKellermann et all (Il977l l used very long baseline interferom¬ 
etry (VLBI) to observe Sgr A* at 7.8 GHz, and they detected 
a nearby secondary transient source which has not been ex¬ 


plaine d, and could in fact be S2 or another star. Further¬ 
more, |Herrnstein]eF^i] ll2004h monitored the flux density of 
Sgr A* using the Very La rge Array. Their results are consis¬ 
tent with our calculation. iMacauart fc Bower! ll2006ll looked 
at the long-term variability of Sgr A* and detected a flux 
of ~ 10®*^ erg s“^ around the 1987 pericenter passage of S2. 
This indicates that the mass loss rate of S2 is ~ 10“® M© 
yr“^. However, these data points are sparse and hence do 
not provide a tight constraint on the peak flux during the 
future pericenter passage. 

While we expect the synchrotron emission to be de¬ 
tectable, the same is not necessarily true for the standoff 
radius. In Figure [3] we plot the standoff radius versus or¬ 
bital radius. The standoff radius near the pericenter may 
not be resolvable, although at larger distances from Sgr A* 
the likelihood for resolving the source increases. It is also 
worth noting that we calculated the thermal free-free emis¬ 
sion around Sgr A* and concluded that it is not detectable. 


4 CONCLUSIONS AND FUTURE 

OBSERVATIONS 

The innermost stars orbiting Sgr A* (such as the S-stars) 
produce a bow shock. Although this shock will likely not be 
detectable via thermal emission, we have shown that stars 
such as S2 may be detectable via synchrotron emission. Fig¬ 
ure [T] shows that the contrast of synchrotron flux relative to 
Sgr A* is maximized around the 1.4 GHz band. However, 
due to scattering we will need to use a higher frequency, 
around 10 GHz, if we are to resolve our star. If a star such 
as S2 emits strong synchrotron emission, the combined sig¬ 
nal from the star and Sgr A* may exceed the quiescent ra¬ 
dio emission from the black hole. As was the case with G2 
(see next paragraph) it is not clear that any additional syn¬ 
chrotron emission will be observable. Arguably, it is best 
to resolve the synchrotron emission from any star orbiting 
Sgr A*, such as S2. At apoapse S2 is ~ 0.23” from Sgr A* 
while at periapse it is only 0.015” away. Thus, to detect 
a star such as S2 requires both good sensitivity and reso¬ 
lution. VLBI can provi de submilliarcsec observations of Sgr 
A* (e.g. ILu et allboilll and thus the required precision. S2 
is an ideal test case since the star will reach periapse in 2017 
or early 2018. 

Recently, gas cloud G2 was observed orbiting Sgr A*, 
and was approximately 3100 Schwarzschild radii from the 
black hole at pericenter. It was predicted that the bow 
shock from this encounter would displace the quiescent 
radio emission of Sgr A* by ~ 33 mas llNaravan et al.l 
I2OI2I : ISadowski et al.l|2013h . However, observations across 
the spectrum showed no apparent variabil i ty during the pe- 
riastr on passage of G2 (iBower et al.ll2of5l : IValencia-S et al] 
I2OI5II . It is unclear why a bow shock was not detected. 
One possibility is that at the cent er of G2 is low-mass star 
with wind velocity ^ 100 km s“^ (iGrumlev &: Kumarj|2013l : 
IScoville fc BurkeTOl2013lf . Synchrotron emission from such a 
small Uu, would be extremely difficult to detect. Star S2 is 
massive, and thus the winds are likely an order of magni¬ 
tude larger. If synchrotron emission is not detected for S2 
it may be that our value for is too large. It may also be 
that the wind speeds from S2 are lower than expected, or 
that the number density around Sgr A* is significantly less 
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Figure 1. Non-therm al synchrotron power and flux compared with emission from Sgr A* (data for Sgr A* was obtained from 
lYuan Naravanll20f4h . The left panel shows the dependence of the synchrotron emission on wind mass loss rate. The star’s veloc¬ 
ity and wind velocity are both fixed at 1000 km s“^. Mass loss, M, was computed with values of IO^^Mq yr“^ (blue line), 10 “®Mq 
yr“^ (green line), and 10 “®Mq yr“^ (red line). In the right panel we kept the mass loss constant at 10 ~®Mq yr“^ and used wind 
velocities of lOOO km s“^(blue line), 2000 km s“^(green line), and 4000 s“^(red line). 




Figure 3. Left: The standoff radius for S2 (Ro) versus the orbital radius (R* in units of 10~® pc). We see that the value is always 
less than unity, as expected. Right: angular diameter of the standoff radius (6o) versus the angular diameter of the orbital radius (6*). 
Around pericenter it may be difficult to resolve the finite size of the standoff radius. However, S2 should still be detectable via synchrotron 
emission. 


than what is predicted. Arguably, the most likely reason why 
synchrotron emission from S2 may not be detected is simply 
that the shock is not as strong as assumed. Our results are 
valid so long as the shock is strong, but the emission could 
be much fainter if the strength of the shock is lessened. Since 
the shock will be strongest when S2 passes periapse, this will 
be the ideal time to monitor the star for synchrotron emis¬ 


sion. Even a null detection wi ll help place constraints on the 
environment around Sgr A* jGiannios fc Sironill2013l) . 

In addition to radio observations at around the 10 GHz 
band. Figure [1] shows that at infrared wavelengths, in par¬ 
ticular around 10^“^ Hz, Sgr A* is approximately the same 
luminosity as a close-by star such as S2. However, S2 itself 
emits in the infrared, therefore we will not be able to distin¬ 
guish between the non-thermal and thermal emission. We 










Detecting Stars at the Galactic Centre via Synchrotron Emission 5 


will need to compare the measured flux with that of the S2. 
An instrument with enough sensitivity, such as the James 
Webb Space Telescope (JWST), should be able to detect the 
synchrotron emission in the infrared. JWST is scheduled to 
launch in October of 2018. While S2 will have passed peri- 
center, if our predictions are correct the synchrotron emis¬ 
sion should still be detectable. Furthermore, it may be pos¬ 
sible to measure the increase in total infrared emission from 
S2 between apocenter and pericenter owing to the increased 
synchrotron emission. While S2 is the litmus test for detect¬ 
ing stars at the Galactic Centre via synchrotron emission, 
we hope that ultimately we can use this technique to detect 
and monitor stars that have thus far not been observable. 
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